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//////////////////////////////////////////////////////////////////////////////// 



// Include Directives 
/ /==================== 

# include 
# include 
# include 
# include 



PRECOMP 

"Analysis/Color/CommonColor .h" 
<stdlib.h> 
<math.h> 



//////////////////////////////////////////////////////////////////////////////// 
// 

// ColorMatch 
// 

// Description: 

// Match the color feature between the input color image and input color feature 

array 

// Feature: 

// the estimated color spectrum 

// Parameters: 

// pSrcRef - source image 

// pRoi - ROI 

// pColSpectrumRef - input spectrum feature 

// pTolerance - match threshold 

// pSatThreshold - threshold for distinguishing two color with the same hue value 

// MatchFlag - Matching indicator 

// MatchScore - Match score 

// Return Value: 

// Error code 

// 

. //////////////////////////////////////////////////////////////////////////////// 



// algorithms for color feature extraction 

enum ColorFeatureSet { 
kColSpectrum = 0, 
kColSignature, 
kColHis Intersection 
}; 



enum ColorComplexity { 
kComplexityLow = 0, 
kComplexityMedium, 
kComplexityHigh 
} ; 



enum ColorBinNumber { 

kColorBinNbLow = 16, 
kColorBinNbMedium=30 , 
kColorBinNbHigh=58 
} ; 
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idefine CoiorSatSpread 0.2 
#define ColorSatRemain 0.8 
idefine BWHueThreshold 10 
#define BlkLgtThreshold 5 
idefine BWLgtMiddle 128 

#define OffsetCclor 18 // the center of the hue value of red color 

GRLIBError ColorMatch (const GRImage* pSrcIMRef , ArraylD* pColSpectrumRef , const ROI* pRoik^ 
, double* ColTolerance, int pSatThreshold, ArraylD* pMatchFlagRef , ArraylD* ^ 
pMatchScoreRef ) { 

GRImage* IMasklmagePtr = NULL; 
MaskPart IMaskPart; 
const Pix8* IMaskPixPtr; 
const Pix8* IMaskLinePtr; 
int lpixelsOutsideMask=0; 

const GRImage* ISrcIMPtr; 

ConstPixelPtr ISrcPixPtr; 

const PixRGB *lSrcPixRGBPtr , *lBuf SrcRGBPtr ; 

const PixHSL *lSrcPixHSLPtr , *lBuf SrcHSLPtr ; 

double *lColImSpecAlDWPtr, *lColTempSpecAlDWPtr ; 

long *lMatchFlagPtr, *lMatchScorePtr; 

// Coord ICoord; 

int i, jj, ISizeXLW, ISizeYLW, ISrcLineWidth, IColLW, ILineLW ; 

double llmageAreaDW, IHueStepDW; 

int IColSpectrumDimLW, INbColFeature, INbHueBin, IHuelndex, >^ 

INiamColorEnt ry=0 ; 

long IRHHVal, IGSSVal, IBLVVal; 

long Ihue, Isat, llgt, IBlkThreshold; 

Byte Ireplace = 191, offset=0; 

long Icoring = 0; 
double h,s,l; 

long IMskLineWidth, INbROIContour, IROIContour Index; 

double IColDif f , IColImSpreaded [ 60] , IColTempSpreaded [ 60] , lColImSpreaded2 [ 60] , \l 
lColTempSpreaded2 [60] ; 

GRLIBError error ^ ERR_SUCCESS; 

// lookup table for function f (x) =exp (-0 . 025*x) , x=a - 200 
static double expLookup[] = { 
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2938, 
2288, 
1782, 

1 -\ n, 
LOOO , 

1081, 
0842, 
0655, 
0510, 
0398, 
0310, 
0241, 
0188, 
0146, 
0114, 
0089, 
0069, 
0054 

); 

// get the pointer to the input image structure 
ISrcIMPtr = pSrcIMRef ; 
ISrcLineWidth = lSrcIMPtr->lineWidth; 
if ( pRoi->nuraContours >0 ) 

lNbROIContour= ( long) pRoi->nuinContours; // ROI connected 
else 

lNbROIContour=l ; // ROI not connected, use the entire image 
// Set the size of match flag array 

if (error = ResizeArraylD (pMatchFlagRef , sizeof ( long) , INbROIContour) ) 
{ 

goto END; 

) 

// Get the match flag array pointer 

if (error = GetArraylDPtr (pMatchFlagRef , &lMatchFlagPtr) ) 
{ 

goto END; 
} 
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// Set the size of match Score array 

if (error = ResizeArraylD (pMatchScoreRef , sizeof (long) , INbROIContour) ) 
{ 

goto END; 

} 

// Get the match score array pointer 

if (error = GetArraylDPtr (pMatchScoreRef , &lMatchScorePtr) ) 
{ 

goto END; 
} 

/*Initilization */ 

for (i=0; i<lNbROIContour ; i++) { 

IMatchScorePtr [i] =0; 
IMatchFlagPtr [ i ] =0 ; 

) 

// Get the color array pointer 

if (error = GetArraylDPtr (pColSpectrumRef, SlColTempSpecAlDWPtr ) ) { 
goto END; 

} 

// Get the size of the input color spec array 

if( error = GetArraylDSize (pColSpectrumRef , &lColSpectrumDimLW) ) { 
goto END; 

} 

if ( (lColImSpecAlDWPtr= (double *) malloc (IColSpectrumDimLW * sizeof (double) ) ) == 
NULL ) { 

er ror=ERR_OUTOFMEMORY ; 
goto END; 

} 

// compute the parameters for the color table 
INbColFeature =lColSpectrumDimLW; 
lNbHueBin={lNbColFeature -2) /2; 
lHueStepDW=255 . 0/ (double) INbHueBin; 

// moving filtering the template's color spectrum 
switch (IColSpectrumDimLW) { 
case kColorBinNbMedium: 
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for (i=2; i<lNbColFeature-4 ; i+=2) { 

lColTempSpreaded[i]= 0.8* IColTempSpecAlDWPtr [i] +0.1* ^£ 
IColTempSpecAlDWPtr [i-2) + 0.1* IColTempSpecAlDWPtr [i+2 ] ; 

lColTeinpSpreaded[i+l]= 0 . 8*lColTempSpecAlDWPtr [i+1 ] + 0.1* ^ 
IColTempSpecAlDWPtr [i-1] + 0.1* IColTempSpecAlDWPtr [i+3 ] ; 

} 

lColTempSpreaded[0]= 0.8* IColTempSpecAlDWPtr [0] +0.1* ^ 
IColTempSpecAlDWPtr [lNbColFeature-4] + 0.1* IColTempSpecAlDWPtr [2 ] ; 

lColTempSpreaded[l]= 0 . 8*lColTempSpecAlDWPtr [ 1 ] +0.1* \£ 
lColTempSpecAlDWPtr[lNbColFeature-3] + 0.1* IColTempSpecAlDWPtr [3] ; 

lColTempSpreaded[lNbColFeature-4]= 0.8* IColTempSpecAlDWPtr ^ 
[lNbColFeature-4] + 0.1* IColTempSpecAlDWPtr [lNbColFeature-6] +0.1* m 
IColTempSpecAlDWPtr [0] ; 

lColTempSpreaded[lNbColFeature-3] = 0 . 8*lColTempSpecAlDWPtr z 
[lNbColFeature-3] +0.1* IColTempSpecAlDWPtr [lNbColFeature-5] +0.1* 
IColTempSpecAlDWPtr [1] ; 

breaks- 
case kColorBinNbHigh: 

for (i-4; i<lNbColFeature-6; i+=2) { 

lColTempSpreaded[i]= 0.7* IColTempSpecAlDWPtr [i ] + 0.1* 
IColTempSpecAlDWPtr [i-2] + 0.1* IColTempSpecAlDWPtr [i+2 ] 

+ 0.05* IColTempSpecAlDWPtr [i-4] + 0.05* >l 

IColTempSpecAlDWPtr [i+4 ] ; 

lColTempSpreaded[i+l]= 0 . 7*lColTempSpecAlDWPtr [i+1 ] + 0.1* )l 
IColTempSpecAlDWPtr [i-1] +0.1* IColTempSpecAlDWPtr [i+3 ] 

+ 0.05* IColTempSpecAlDWPtr [i-3] + 0.05* M 

IColTempSpecAlDWPtr [i+5] ; 
} 

lColTempSpreaded[0]= 0.7* IColTempSpecAlDWPtr [0 ] + 0.1* ^ 
IColTempSpecAlDWPtr [lNbColFeature-4] +0.1* IColTempSpecAlDWPtr [2 ] 

+ 0.05* IColTempSpecAlDWPtr [lNbColFeature-6] + 0. \i 

05* IColTempSpecAlDWPtr [4] ; 

lColTempSpreaded[l]= 0 . 7*lColTempSpecAlDWPtr [ 1 } +0.1* ^ 
IColTempSpecAlDWPtr [lNbColFeature-3] +0.1* IColTempSpecAlDWPtr [3] 

+ 0.05* IColTempSpecAlDWPtr [lNbColFeature-5] + k^' 

0.05* IColTempSpecAlDWPtr [5] ; 

lColTempSpreaded[2]= 0.7* IColTempSpecAlDWPtr [2 ] +0.1* ^ 
IColTempSpecAlDWPtr [0] +0.1* IColTempSpecAlDWPtr [4 ] 

+ 0.05* IColTempSpecAlDWPtr [lNbColFeature-4]+ 0. \t 

05* IColTempSpecAlDWPtr [6] ; 

lColTempSpreaded[3]= 0 . 7*lColTempSpecAlDWPtr [3] +0.1* ^ 
IColTempSpecAlDWPtr [1] +0.1* IColTempSpecAlDWPtr [5] 

+ 0.05* IColTempSpecAlDWPtr [lNbColFeature-3]+i^ 

0.05* IColTempSpecAlDWPtr [7] ; 

lColTempSpreaded[lNbColFeature-4]= 0.7* IColTempSpecAlDWPtr ^ 
[lNbColFeature-4] + 0.1* IColTempSpecAlDWPtr [lNbColFeature-6] +0.1* \t 
IColTempSpecAlDWPtr [0] 
+ 0.05* IColTempSpecAlDWPtr [lNbColFeature-8]+ 0.05* IColTempSpecAlDWPtr [2 ] ; 
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IColTempSpreaded [ lNbCoiFeature-3 ] = 0.7* IColTempSpecAl DWPt r wr 
[lNbColFeature-3] + 0.1* IColTempSpecAlDWPtr [lNbColFeature-5] + 0.1* »^ 
IColTempSpecAlDWPtr [1] 
+ 0.05* IColTempSpecAlDWPtr [lNbColFeature-7]+ 0.05* IColTempSpecAlDWPtr [3 ] ; 

IColTempSpreaded [lNbColFeature-6]= 0.7* IColTempSpecAlDWPtr 
[lNbCoiFeature-6] + 0.1* IColTempSpecAlDWPtr [lNbColFeature-8 ] + 0.1* y£ 
IColTempSpecAlDWPtr [lNbColFeature-4] 
+ 0.05* IColTempSpecAlDWPtr [lNbColFeature-10] + 0.05* IColTempSpecAlDWPtr [ 0] ; 

IColTempSpreaded [lNbColFeature-5]= 0 . 7*lColTempSpecAlDWPtr si 
tlNbColFeature-5] +0.1* IColTempSpecAlDWPtr [lNbColFeature-7 ] +0.1* s£ 
IColTempSpecAlDWPtr tlNbColFeature-3] 
+ 0.05* IColTempSpecAlDWPtr [lNbColFeature-9]+ 0.05* IColTempSpecAlDWPtr [ 1 ] ; 
break; 

case JcColorBinNbLow : 
default : 

for (i=0; i<lNbColFeature-2 ; i+=2) { 

IColTempSpreaded [ i ] = IColTempSpecAlDWPtr [ i ] ; 
IColTempSpreaded [i+l]= IColTempSpecAlDWPtr [i+1] ; 

} 



break; 



} 



// fuzzy a little bit for saturation too 

for (i^O; i<lNbColFeature-2 ; i+=2) { 

lColTempSpreaded2 [i]= ColorSatRemain*lColTempSpreaded [i] + ColorSatSpread * ^ 
IColTempSpreaded [i+1] ; 

lColTempSpreaded2 [i+1] = ColorSatRemain* IColTempSpreaded [i+1 ] + ColorSatSpreadi^ 
* IColTempSpreaded [ i ] ; 

} 

// no fuzzy for black and white 

lColTempSpreaded2 [ lNbColFeature-2 ] =lColTempSpecAl DWPtr [ lNbColFeature-2 ] ; 
lColTempSpreaded2 [lNbColFeature-1 ] =lColTempSpecAl DWPtr [lNbColFeature-1] ; 

// if ROI input connected, calculate the color feature based on the ROI image 
if ( pRoi->numContours>0 ) { 

// calculate the color match for each subROI 

IROIContourlndex^O; • 

j j =lNbROICont our ; 
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while(jj — ) { 

/*Initilization */ 
for (i=0; KlNbColFeature; { 
IColImSpecAlDWPtr [i] =0; 

} 

// create a temporary mask image 

if ( error= New Image ( &lMaskImagePtr , IMAGE_U8, 4, 4, 0) ) { 
goto END; 

} 

// ROIContourToMask (const ROI* _roi, int _contourNdx, int _fillValue, GRImage* 
_mask Image) ; 

if (error =ROIContourToMask (pRoi, IROIContourlndex, 1, IMasklmagePtr) ) { 
goto END; 

}; 

IMskLineWidth = lMaskImagePtr->lineWidth; 
CalcMaskPart (pSrcIMRef , IMasklmagePtr, &lMaskPart) ; 

GetConstlmagePixel (pSrcIMRef , IMaskPart . startlmageX, IMaskPart . startlmageY, & 
ISrcPixPtr) ; 

GetConstlmagePixel (IMasklmagePtr, IMaskPart . startMaskX, IMaskPart . startMaskY, 
(ConstPixelPtr* ) &lMaskLinePtr ) ; 

ISizeXLW = IMaskPart .sizeX; 
ISizeYLW = IMaskPart . sizeY; 

llmageAreaDW =0; /* init the Area */ 



// control negative sizes 

if (ISizeXLW < 0) 

ISizeXLW - 0; 
if (ISizeYLW < 0) 

ISizeYLW = 0; 



switch (lSrcIMPtr->imageType) { 
case IMAGE_RGB32: 

lSrcPixRGBPtr=lSrcPixPtr : PixRGB_Ptr; 
break; 
case IMAGE_HSL32: 

lSrcPixHSLPtr=lSrcPixPtr . PixHSL Ptr ; 
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break; 

default : 

error = ERR_BADIMAGETYPE; 
goto END; 

} 



// main switch 

switch {lSrcIMPtr->imageType) { 

case IMAGE_RGB32: 

for (ILineLW = 0; ILineLW < ISizeYLW; lLineLW++)' { 

lBufSrcRGBPtr=lSrcPixRGBPtr; 
IMaskPixPtr = IMaskLinePtr ; 

for {IColLW = 0; IColLW < ISizeXLW; IC0ILW++) { 
if (*lMaskPixPtr++) { 



// initializatin 

lNumColorEntry=0 ; 

// convert the RGB to HSL that corresponds to our board 

IRHHVal = lBufSrcRGBPtr->Byte_RGB.R; 
IGSSVal = lBufSrcRGBPtr->Byte_RGB.G; 
IBLVVal = lBufSrcRGBPtr->Byte_RGB.B; 

CCONV_HSI (IRHHVal, IGSSVal , IBLVVal , Icoring, lreplace,l, s,h) 

Ihue = (Byte) (h + 0.5 + offset) % 256; 
lsat= (Byte) (s+0.5); 
llgt = (Byte) (1+0.5) ; 
// clustering the color feature 
if( lsat< BWHueThreshold ) { 

lNumColorEntry=(llgt < BWLgtMiddle ) ? (lNbColFeature-2) : »^ 

(lNbColFeature-1) ; 

} 



else { 

// Threshold Black color using Nobelist*s formular 

// lBlkThreshold= (long) ((128 - BlkLgtThreshold) *exp (-0 . 
025*(lsat - BWHueThreshold)) + BlkLgtThreshold ); 

ifdsat < 200 ) 

lBlkThreshold= (long) ((128 - BlkLgtThreshold) * 
expLookup[lsat - BWHueThreshold] + BlkLgtThreshold ); 

else 

lBlkThreshold=BlkLgtThreshold; 
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if {llgt < IBlJcThreshold ) { 

lNuniColorEntry= ( lNbColFeature-2 ) ; 

} 

else { 

// offset the red color to the center of the first 

Ihue = ( Ihue + OffsetCoior) % 255; 
// clustering non-black-white color feature 
lHueIndex= (double) Ihue / IHueStepDW; 
lNumColorEntry= (Isat < pSatThreshold) ? (IHuelndex* 

} 

} 

IColImSpecAlDWPtr [INumColorEntry] += 1.0; 

lImageAreaDW++; 
} // end of if (*lMaskPixPtr++) 

IBuf SrcRGBPtr++ ; 

} 

ISrcPixRGBPtr +=lSrcLineWidth; 
IMaskLinePtr += IMskLineWidth; 

} 

// calculate the discrete color spectrum 
// just in case, llmageAreaDW should never be 0 
if (lIinageAreaDW==0 ) 

for( i=0; i< INbColFeature; i++ ) 
IColImSpecAlDWPtr [i] =0; 

else 

for( i=0; i< INbColFeature; i++ ) 

IColImSpecAlDWPtr [i] /= llmageAreaDW; 

break; 
case IMAGE_HSL32: 

for (ILineLW = 0; ILineLW < ISizeYLW; lLineLW++) { 



bin 



2) : (lHueIndex*2+l) ; 
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lBufSrcHSLPtr=lSrcPixHSLPtr; 
IMaskPixPtr = IMaskLinePtr ; 

for (IColLW = 0; IColLW < ISizeXLW; IC0ILW++) { 

if {*lMaskPixPtr++) { 

// initializatin 

lhue=lBuf SrcHSLPtr->Byte_HSL . H; 
lsat=lBufSrcHSLPtr->Byte_HSL.S; 
llgt=lBuf SrcHSLPtr->Byte_HSL . L; 

lN\ainColorEntry=0 ; 

if( lsat< BWHueThreshold ) { 

lNuitiColorEntry={llgt < BWLgtMiddle ) ? (lNbColFeature-2) : 

(lNbColFeature-1) ; 

} 

else { 

// Threshold Black color using Nobelist's formular 

// lBlkThreshold= (long) ((128 - BlkLgtThreshold) *exp (-0 . 025* 
(Isat - BWHueThreshold)) + BlkLgtThreshold }; 

ifdsat < 200 ) 

lBlkThreshold= (long) ((128 - BlkLgtThreshold) * 
expLookup[lsat - BWHueThreshold] + BlkLgtThreshold ); 

else 

lBlkThreshold=BlkLgtThreshold; 
if (llgt < IBlkThreshold ) { 

lNumColorEntry= (lNbColFeature-2 ) ; 

} 

else { 

// offset the red color to the center of the first bin 
Ihue = ( Ihue + OffsetColor) % 255; 
// clustering non-black-white color feature 
lHueIndex= (double) Ihue / IHueStepDW; 

lNumColorEntry= (Isat < pSatThreshold) ? (lHueIndex*2) : 

{lHueIndex*2+l) ; 

} 

} 



IColImSpecAlDWPtr [INumColorEntry] += 1.0; 
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lImageAreaDW++; 

} // end of if {*lMaskPixPtr++) 
lBufSrcHSLPtr++; 



lSrcPixHSLPtr+=lSrcLineWidth; 
IMaskLinePtr += IMskLineWidth; 

} 

// calculate the discrete color spectrum 
// just in case, llmageAreaDW should never be 0 
if (lImageAreaDW==0 ) 

for{ i=0; i< INbColFeature; i++ ) 
IColImSpecAlDWPtr [i] =0; 

else 

for( i=0; i< INbColFeature; i++ ) 

IColImSpecAlDWPtr [i] /= llmageAreaDW; 

break; 

} // end of switch (lSrcIMPtr->imageType) 

// moving filtering input image's color spectrum 

switch (IColSpectrumDimLW) { 

case kColorBinNbMedium: 

for {i=2; i<lNbColFeature-4 ; i+=2) { 

lColImSpreaded[i] = 0.8* IColImSpecAlDWPtr [i] + 0.1* 
IColImSpecAlDWPtr [i-2] +0.1* IColImSpecAlDWPtr [i+2] ; 

lColImSpreaded[i+l] = 0.8* IColImSpecAlDWPtr [i+1 ] + 0.1* M 
IColImSpecAlDWPtr [i-1] +0.1* IColImSpecAlDWPtr [i+3] ; 

} 

lColImSpreaded[0] = 0.8* IColImSpecAlDWPtr [ 0] +0.1* 
IColImSpecAlDWPtr [lNbColFeature-4] + 0.1* IColImSpecAlDWPtr [2 ] ; 

lColImSpreaded[l] = 0.8* IColImSpecAlDWPtr [ 1] + 0.1* M 
IColImSpecAlDWPtr [lNbColFeature-3] +0.1* IColImSpecAlDWPtr [3] ; 

lColImSpreadedtlNbColFeature-4] =0.8* IColImSpecAlDWPtr 
[lNbColFeature-4] + 0.1* IColImSpecAlDWPtr [lNbColFeature-6] + 0.1* IColImSpecAlDWPtr 
[0]; 

lColImSpreaded[lNbColFeature-3] =0.8* IColImSpecAlDWPtr \£ 
[lNbColFeature-3] +0.1* IColImSpecAlDWPtr [lNbColFeature-5] +0.1* IColImSpecAlDWPtr 
[1]; 
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break; 

case kColorBinNbHigh : 

for (i=4; i<lNbColFeature-6; i+=2) { 

lColImSpreaded[i] =0.7* IColImSpecAlDWPtr [ i] +0.1* 
IColImSpecAlDWPtr [i-2] +0.1* IColImSpecAlDWPtr [i+2 ] 

+ 0.05* IColImSpecAlDWPtr [i-4] +0.05* M 

IColImSpecAlDWPtr [i+4] ; 

lColImSpreaded[i+l] = 0.7* IColImSpecAlDWPtr [i+1] + 0.1* y£ 
IColImSpecAlDWPtr [i-1] +0.1* IColImSpecAlDWPtr [i+3] 

+ 0.05* IColImSpecAlDWPtr [i-3] + 0.05* \£ 

IColImSpecAlDWPtr [i+5] ; 



) 



lColImSpreaded[0] =0.7* IColImSpecAlDWPtr [0] +0.1* )t 
IColImSpecAlDWPtr [lNbColFeature-4] +0.1* IColImSpecAlDWPtr [2 ] 

+ 0.05* IColImSpecAlDWPtr [lNbColFeature-6]+ O.OSi^ 

* IColImSpecAlDWPtr [4] ; 

lColImSpreaded[l] = 0.7* IColImSpecAlDWPtr [ 1 ] +0.1* ^ 
IColImSpecAlDWPtr [lNbColFeature-3] +0.1* IColImSpecAlDWPtr [3] 

+ 0.05* IColImSpecAlDWPtr [lNbColFeature-5]+ 0.05^ 

* IColImSpecAlDWPtr [5] ; 

lColImSpreaded[2] =0.7* IColImSpecAlDWPtr [2 ] +0.1* M 
IColImSpecAlDWPtr [0] +0.1* IColImSpecAlDWPtr [4 ] 

+ 0.05* IColImSpecAlDWPtr [lNbColFeature-4)+ 0.05»^ 

* IColImSpecAlDWPtr [6] ; 



lColImSpreaded[3] = 0.7* IColImSpecAlDWPtr [3 ] + 0.1* \£ 
IColImSpecAlDWPtr [1] +0.1* IColImSpecAlDWPtr [5] 

+ 0.05* IColImSpecAlDWPtr [lNbColFeature-3]+ 0.05»*r 

* IColImSpecAlDWPtr [7] ; 



lColImSpreaded[lNbColFeature-4] =0.7* IColImSpecAlDWPtr >l 
[lNbColFeature-4] + 0.1* IColImSpecAlDWPtr [lNbColFeature-6] +0.1* IColImSpecAlDWPtr 
[0] 

+ 0.05* IColImSpecAlDWPtr ^ 

[lNbColFeature-8)+ 0.05* IColImSpecAlDWPtr [2 ] ; 

lColImSpreaded[lNbColFeature-3] =0.7* IColImSpecAlDWPtr ^ 
[lNbColFeature-3] +0.1* IColImSpecAlDWPtr [lNbColFeature-5] +0.1* IColImSpecAlDWPtr M 
[1] 

+ 0.05* IColImSpecAlDWPtr [lNbColFeature-7]+ 0.05* IColImSpecAlDWPtr [3] ; 



lColImSpreaded[lNbColFeature-6] = 0.7* IColImSpecAlDWPtr ^ 

[lNbColFeature-6] + 0.1* IColImSpecAlDWPtr [lNbColFeature-8] +0.1* IColImSpecAlDWPtr >l 
tlNbColFeature-4) 

+ 0.05* IColImSpecAlDWPtr M 

[lNbColFeature-10]+ 0.05* IColImSpecAlDWPtr [0] ; 

lColImSpreaded[lNbColFeature-5] =0.7* IColImSpecAlDWPtr \i 

[lNbColFeature-5] + 0.1* IColImSpecAlDWPtr [lNbColFeature-7 ] +0.1* IColImSpecAlDWPtr ^ 
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[lNbColFeature-3] 

+ 0.05* IColImSpecAlDWPtr [lNbColFeature-9]+ 0.05* IColImSpecAlDWPtr [1] ; 
break; 

case kColorBinNbLow : 
default: 

for (i=0; i<lNbColFeature-2 ; i+=2) { 

lColImSpreaded[i] = IColImSpecAlDWPtr [i] ; 
lColImSpreaded[i+l] = IColImSpecAlDWPtr [i+1 ] ; 

} 

break; 



// fuzzy a little bit for saturation too 
for (i=0; i<lNbColFeature-2 ; i+=2) { 

lColImSpreaded2 [i] =ColorSatRemain*lColImSpreaded [i] + ColorSatSpread* 
lColImSpreaded[i+l] ; 

lColImSpreaded2 [i+1] =ColorSatRemain*lColImSpreaded [i+1 ] + ColorSatSpread * ^ 
lColImSpreaded[i] ; 

} 

// no fuzzy for black and white 
lColImSpreaded2 [lNbColFeature-2 ] 
lColImSpreaded2 [lNbColFeature-1 ] 

// calculate the color distance 
lColDiff=0; 
for (i=0; KlNbColFeature; i++) { 

IColDiff += fabs (lColTempSpreaded2 [i] - lColImSpreaded2 [i] ); 

// histogram instersection if # of bin is large enough 

// IColDiff += (lColTempSpreaded2[i] <= lColImSpreaded2 [i] ) ? »^ 
lColTempSpreaded2 [i] : lColImSpreaded2 [i] ; 



=lColImSpecAlDWPtr [lNbColFeature-2] ; 
= IColImSpecAlDWPtr [lNbColFeature-1] ; 



} 

if (lImageAreaDW==0 ) // this should never happen ! 
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lColDiff=4; 

if ( IColDiff < *ColTolerance *2.0 ) 

IMatchFlagPtr [lROIContourIndex]=l; 
else 

IMatchFlagPtr [lROIContourIndex]=0; 
lMatchScorePtr [IROIContourlndex] = (2 - IColDiff) *500; // get the magic 
score number 

/ /des tory the temporary mask image 
if (IMasklmagePtr ! =NULL ) 

Disposelmage (IMasklmagePtr) ; 
// get next contour index 
IROIContour Index++ ; 

} // end of while (INbROIContour--) 

} // end of ROI connected case 

else { // ROI not connected, take care of the entire image now 
/*Initilization */ 
for {i=0; i<lNbColFeature; i++) { 
IColImSpecAlDWPtr [i]=0; 

} 

// there is only one element in the match flag and match score array 
IROIContour lndex=0 ; 

GetConstlmagePixel (pSrcIMRef , 0,0, ilSrcPixPtr) ; 

ISizeXLW = lSrcIMPtr->xRes; 
ISizeYLW = lSrcIMPtr->yRes; 

llmageAreaDW = (double) ISizeXLW * ISizeYLW; 
switch (lSrcIMPtr->imageType) { 
case IMAGE_RGB32: 

lSrcPixRGBPtr=lSrcPixPtr.PixRGB_Ptr; 
break; 
case IMAGE_HSL32 : 

lSrcPixHSLPtr=lSrcPixPtr.PixHSL Ptr; 
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breaks- 
default: 

error = ERR_BADIMAGETYPE; 
goto END; 

} 

// main switch 

switch (ISrcIMPtr-MmageType) { 
case IMAGE_RGB32: 

for (ILineLW = 0; ILineLW < ISizeYLW; lLineLW++) { 
lBufSrcRGBPtr=lSrcPixRGBPtr; 

for (IColLW = 0; IColLW < ISizeXLW; IC0ILW++) { 

// initializatin 

INiimColorEntry^O; 

// convert the RGB to HSL that corresponds to our board 

IRHHVal = lBufSrcRGBPtr->Byte_RGB.R; 
IGSSVal = lBufSrcRGBPtr->Byte_RGB.G; 
IBLVVal = lBufSrcRGBPtr->Byte_RGB.B; 

CCONV_HSI (IRHHVal, IGSSVal, IBLVVal, Icoring, Ireplace, 1, s, h) 

Ihue = (Byte) (h + 0.5 + offset) % 256; 
lsat= (Byte) (s+0.5); 
llgt = (Byte) (1+0.5) ; 
// clustering the color feature 
if( lsat< BWHueThreshold ) { 

lNuinCoiorEntry=(llgt < BWLgtMiddle ) ? (lNbColFeature-2) : 

} 

else { 

// Threshold Black color using Nobelist's formular 

// lBlkThreshold= (long) ((128 - BlkLgtThreshold) *exp (-0 , »^ 
025*(lsat - BWHueThreshold)) + BlkLgtThreshold ); 

ifdsat < 200 ) 

lBlkThreshold= (long) ((128 - BlkLgtThreshold) * 
expLookup [Isat - BWHueThreshold] + BlkLgtThreshold ); 

else 

lBlkThreshold=BlkLgtThreshold; 
if (llgt < IBlkThreshold ) { 

lNumColorEntry= (lNbColFeature-2 ) ; 



(lNbColFeature-1) 
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} 

else { 

// offset the red color to the center of the first ^ 

bin 

Ihue = ( Ihue + OffsetColor) % 255; 
// clustering non-black-white color feature 
lHueIndex= (double) Ihue / IHueStepDW; 
lNumColorEntry= (Isat < pSatThreshold) ? (IHuelndex* 

2) : (lHueIndex*2+l) ; 

} 

} 

IColImSpecAlDWPtr [INumColorEntry] 1.0/ 
lBufSrcRGBPtr++; 

} 

ISrcPixRGBPtr +=lSrcLineWidth; 

} 

// calculate the discrete color spectrum 
// just in case, llmageAreaDW should never be 0 
if (lImageAreaDW==0 ) 

for( 1=0; i< INbColFeature; i++ ) 
IColImSpecAlDWPtr [i] =0; 

else 

for( i=0; i< INbColFeature; i++ ) 

IColImSpecAlDWPtr [i] /= llmageAreaDW; 

break; 
case IMAGE_HSL32: 

for (iLineLW = 0; iLineLW < ISizeYLW; lLineLW++) { 
lBufSrcHSLPtr=lSrcPixHSLPtr; 

for (IColLW = 0; IColLW < ISizeXLW; IC0ILW++) { 

// initial iza tin 

lhue=lBuf SrcHSLPtr->Byte_HSL . H; 
lsat=lBufSrcHSLPtr->Byte_HSL.S; 
llgt=lBuf SrcHSLPtr->Byte_HSL. L; 

INumColorEntry^O; 
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if( lsat< BWHueThreshold ) { 

lNumColorEntry=(llgt < BWLgtMiddle ) ? {lNbColFeature-2) : >i 

(INbColFeature-l) ; 

} 

else { 

// Threshold Black color using Nobelist's formular 

// lBlkThreshold= (long) ((128 - BlkLgtThreshold) *exp (-0 . 025* »^ 
(Isat - BWHueThreshold)) + BlkLgtThreshold ); 

if (Isat < 200 ) 

lBlkThreshold= (long) ((128 - BlkLgtThreshold) * M 
expLookup [Isat - BWHueThreshold] + BlkLgtThreshold ); 

else 

lBlkThreshold=BlkLgtThreshold; 
if (llgt < IBlkThreshold ) { 

lNumColorEntry= ( lNbColFeature-2 ) ; 

} 

else { 

// offset the red color to the center of the first bin 
Ihue = ( Ihue + OffsetColor) % 255; 
// clustering non-black-white color feature 
lHueIndex= (double) Ihue / IHueStepDW; 

lNuinColorEntry= (Isat < pSatThreshold) ? (lHueIndex*2 ) : \l 

(lHueIndex*2+l) ; 

} 

1 

IColImSpecAlDWPtr [INumColorEntry] += 1.0; 
lBufSrcHSLPtr++; 

} 

lSrcPixHSLPtr+=lSrcLineWidth; 

} 

// calculate the discrete color spectrum 
// just in case, llmageAreaDW should never be 0 

if (lImageAreaDW==0 ) 

for( i=0; i< INbColFeature; i++ ) 

IColImSpecAlDWPtr [i] =0; 

else 
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for{ i=0; i< INbColFeature; i++ ) 

IColImSpecAlDWPtr [i] /= llmageAreaDW; 

break; 

} // end of switch (lSrcIMPtr->imageType) 

// moving filtering input image's color spectrum 

switch (IColSpectrumDimLW) { 

case kColorBinNbMedium: 

for {i=2; i<lNbColFeature-4 ; i+=2) { 

lColImSpreaded[i] = 0.8* IColImSpecAlDWPtr [i] + 0.1* 
IColImSpecAlDWPtr [i-2] +0.1* IColImSpecAlDWPtr [i+2 ] ; 

lColImSpreaded[i-hl] = 0.8* IColImSpecAlDWPtr [i+1] + 0.1* 
IColImSpecAlDWPtr [i-1] +0.1* IColImSpecAlDWPtr [i+3] ; 

} 

lColImSpreaded[0] = 0.8* IColImSpecAlDWPtr [0] + 0.1* 
IColImSpecAlDWPtr [lNbColFeature-4] + 0.1* IColImSpecAlDWPtr [2 ] ; 



lColImSpreaded[l] = 0.8* IColImSpecAlDWPtr [1] + 0.1* wr 
IColImSpecAlDWPtr [lNbColFeature-3] + 0.1* IColImSpecAlDWPtr [3] ; 

lColImSpreaded[lNbColFeature-4] = 0.8* IColImSpecAlDWPtr 
[lNbColFeature-4] +0.1* IColImSpecAlDWPtr [lNbColFeature-6] +0.1* IColImSpecAlDWPtr \t 
[0]; 

lColImSpreaded[lNbColFeature-3] = 0.8* IColImSpecAlDWPtr \i 
[lNbColFeature-3] +0.1* IColImSpecAlDWPtr [lNbColFeature-5] +0.1* IColImSpecAlDWPtr yl 
[1]; 



breaks- 
case kColorBinNbHigh : 

for (i=4; i<lNbColFeature-6; i+=2) { 

lColImSpreaded[i] =0.7* IColImSpecAlDWPtr [i] +0.1* 
IColImSpecAlDWPtr [i-2] + 0.1* IColImSpecAlDWPtr [i+2 ] 

+ 0.05* IColImSpecAlDWPtr [i-4] +0.05* 

IColImSpecAlDWPtr [i+4] ; 

lColImSpreaded[i+l] =0.7* IColImSpecAlDWPtr [i+1 ] +0.1* M 
IColImSpecAlDWPtr [i-1] + 0.1* IColImSpecAlDWPtr [i+3] 

+ 0.05* IColImSpecAlDWPtr [i-3] +0.05* ^t 

IColImSpecAlDWPtr [i+5] ; 



} 



lColImSpreaded[0] = 0.7* IColImSpecAlDWPtr [ 0] + 0.1* 
IColImSpecAlDWPtr [lNbColFeature-4] +0.1* IColImSpecAlDWPtr [2 ] 
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+ 0.05* lColImSpecAlDWPtr[lNbColFeature-6]+ O.OSk^ 

* IColImSpecAlDWPtr [4] ; 

lColImSpreaded[l] =0.7* IColImSpecAlDWPtr [ 1] +0.1* 
IColImSpecAlDWPtr [lNbColFeature-3] + 0.1* IColImSpecAlDWPtr [3 ) 

+ 0.05* IColImSpecAlDWPtr [lNbColFeature-5]+ 0.05*^ 

* IColImSpecAlDWPtr [5] ; 

lColImSpreaded[2] = 0.7* IColImSpecAlDWPtr [2] + 0.1* \£ 
IColImSpecAlDWPtr [0] + 0.1* IColImSpecAlDWPtr [4 ] 

+ 0.05* IColImSpecAlDWPtr [lNbColFeature-4]+ 0.05»^' 

* IColImSpecAlDWPtr [6] ; 



lColImSpreaded[3] = 0.7* IColImSpecAlDWPtr [ 3] + 0.1* \^ 
IColImSpecAlDWPtr [1] + 0.1* IColImSpecAlDWPtr [5] 

+ 0.05* IColImSpecAlDWPtr [lNbColFeature-3]+ O.OSi^^ 

* IColImSpecAlDWPtr [7] ; 



lColImSpreaded[lNbColFeature-4] =0.7* IColImSpecAlDWPtr 
[lNbColFeature-4] + 0.1* IColImSpecAlDWPtr [INbColFeature- 6] + 0.1* IColImSpecAlDWPtr 
[0] 

+ 0.05* IColImSpecAlDWPtr wr 

[lNbColFeature-8]+ 0.05* IColImSpecAlDWPtr [2 ] ; 

lColImSpreaded[lNbColFeature-3] =0.7* IColImSpecAlDWPtr 
[lNbColFeature-3] + 0.1* IColImSpecAlDWPtr [lNbColFeature-5] +0.1* IColImSpecAlDWPtr \l 
[1] 

+ 0.05* IColImSpecAlDWPtr [lNbColFeature-7]+ 0.05* IColImSpecAlDWPtr [3 ] ; 

lColImSpreaded[lNbColFeature-6] =0.7* IColImSpecAlDWPtr 
[lNbColFeature-6] +0.1* IColImSpecAlDWPtr [lNbColFeature-8 ] +0.1* IColImSpecAlDWPtr )i 
[lNbColFeature-4] 

+ 0.05* IColImSpecAlDWPtr )^ 

[lNbColFeature-10]+ 0.05* IColImSpecAlDWPtr [0] ; 

lColImSpreaded[lNbColFeature-5] =0.7* IColImSpecAlDWPtr \^ 
[lNbColFeature-5] +0.1* IColImSpecAlDWPtr [lNbColFeature-7 ] +0,1* IColImSpecAlDWPtr )i 
[lNbColFeature-3] 

+ 0.05* IColImSpecAlDWPtr [lNbColFeature-9]+ 0.05* IColImSpecAlDWPtr [ 1] ; 
break; 

case kColorBinNbLow: 
default : 

for (i=0; KlNbColFeature-2 ; i+=2) { 

lColImSpreaded[i] = IColImSpecAlDWPtr [ i] ; 
lColImSpreaded[i+l] = IColImSpecAlDWPtr [i+1] ; 



} 



break; 



} 
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// fuzzy a little bit for saturation too 
for (i=0; i<lNbColFeature-2 ; i+=2) { 

iColIruSpreaded2 [i] =ColorSatReinain*lColImSpreaded [i ] + ColorSatSpread* »^ 
lColImSpreaded[i+l] ; 

lColIinSpreaded2 [i+1] =ColorSatRemain*lColImSpreaded [i+1] + ColorSatSpread * m 
lColImSpreaded[i] ; 

} 

// no fuzzy for black and white 

lColImSpreaded2 [lNbColFeature-2] =lColImSpecAlDWPtr [lNbColFeature-2 ] ; 

lColImSpreaded2 [lNbColFeature-1] = IColImSpecAlDWPtr [ lNbColFeature-1 ] ; 

// calculate the color distance 
lColDiff=0; 
for. (i=0; KlNbColFeature; i++) { 

IColDiff += fabs (lColTeinpSpreaded2 [i] - lColImSpreaded2 [i] ); 

// histogram instersection if # of bin is large enough 

// IColDiff += (lColTempSpreaded2 [i] <= lColImSpreaded2 [i] ) ? 
lColTempSpreaded2 [i] : lColImSpreaded2 [i ] ; 

} 

if (lImageAreaDW==0 ) // this should never happen ! 
IColDiff =4; 

if { IColDiff < *ColTolerance *2.0 ) 
IMatchFlagPtr [IROIContourlndex] =1; 

else 

IMatchFlagPtr [IROIContourlndex] =0; 

IMatchScorePtr [IROIContourlndex] = (2 - IColDiff) *500; // get the magic 
score number 

} // end of if { pRoi->numContours>0) 



END: 

// just do a cleanup 
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if (IColImSpecAlDWPtr ! =NULL) 
free(lColImSpecAlDWPtr) ; 

return error; 

} 



